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Abstract: 

We study two dimensional freely decaying magnetohydrodynamic turbulence. We 
investigate the time evolution of the probability law of the gauge field and the stream 
function. Assuming that this probability law is initially defined by a statistical field theory 
in the basin of attraction of a renormalisation group fixed point, we show that its time 
evolution is generated by renormalisation transformations. In the long time regime, the 
probability law is described by non-unitary conformal field theories. In that case, we prove 
that the kinetic and magnetic energy spectra are proportional. We then construct a family 
of fixed points using the (p, p + 2) non-unitary minimal models of conformal field theories. 
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It has been recently noticed M by Polyakov that field theoretic methods can be used 
to study steady states of two dimensional turbulence. In this approach, Polyakov assumes 
that the velocity field of a turbulent fluid can be described by a primary field deduced from 
a conformal field theory. The same method has been applied^ 2 ' by Ferretti and Yang in the 
case of steady states of magnetohydro dynamic turbulence. Other possible applications to 
magnetohydrodynamic have been also discussed by Coceal and Thomas^. One of the re- 
sults found by Ferretti and Yang is that some steady states do not respect the equipartition 
W between the magnetic and the kinetic energies. It has been argued^ 5 ! by Rahimi-Tabar 
and Rouhani that the equipartition can be satisfied provided that time as well as spa- 
tial scale invariances hold. In this communication, we shall use renormalisation group 
techniques to study the long time behaviour of freely decaying magnetohydrodynamic tur- 
bulence in two dimensions. In particular, we shall see that conformal field theories appear 
as fixed points of the time evolution of the probability law of the stream function and the 
gauge field in the long time regime. At these fixed points, the kinetic and the magnetic 
energy spectra are proportional. Finally, we shall construct a family of fixed points using 
the (p, p + 2) non-unitary minimal models. 

In two dimensions, the velocity field of an incompressible fluid is determined by a 
single pseudo-scalar field: the stream function 

v a = eapdpip. (1) 

The tensor e a p is the antisymmetric tensor e\2 = —e2\ = 1- The vorticity is then given by 

w = e a pd a vp. (2) 

It is a pseudo-scalar too. We are interested in the two dimensional situation where the 
magnetic field lies in the plane of the fluid motion. Moreover we suppose that there is 
no external magnetic field and the magnitude of the initial magnetic field is of the same 
order as the magnitude of the initial velocity field (in well-chosen units). This situation 
can occur for certain astrophysical phenomena where the magnetic Reynolds number is 
large. In the absence of magnetic monopoles, the magnetic field can also be represented 
by 

B a = eapdpA, (3) 
where A is a scalar gauge field. The induced electric current can be deduced 

J = e a pd a Bp. (4) 

These quantities satisfy coupled partial differential equations. We will focus on the free 
decaying case in an infinite plane. We shall suppose that the vorticity field and the electric 
current decay sufficiently fast at infinity. 

The Navier-Stokes equations are effective equations describing the behaviour of a fluid 
in the hydrodynamic approximation. The velocity field is the average of the microscopic 
velocity of molecules over a macroscopic size a, the size of a fluid element. Similarly, the 
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magnetic field is also a macroscopic field. We shall analyse the macroscopic equations for 
the regularised vorticity and gauge fields ^ 



d t uj a + e a pd a ip a dpAij; a = e ct f 3 d a J a d ( 3A a + isAu a 
d t A a + v aa d a A a = 7]J a . 

The viscosity v and the resistivity r\ are responsible for dissipation. As the resistivity and 
the viscosity have the same dimensions (in well-chosen units) , one deduces 

V = Vo — , 

T 

where r]o is of order one. The time r is a characteristic time for the decay of both the 
kinetic and the magnetic energies. 

In field theoretic terms, the size a plays the role of an intrinsic ultraviolet cut-off. 
We shall be interested in the situation where the gauge field and the stream function 
are random fields regularised at short distances by the molecular size a. As no random 
stirring is imposed, we shall suppose that initial conditions are randomly prescribed by a 
probability law for both the initial magnetic and velocity fields. The evolution of a given 
configuration is given by (5). The evolution of the probability law of the gauge field and 
the stream function dV(A a , ip a ) is governed by the Hopf equations obtained considering (5) 
as valid in any equal time correlation functions of the gauge field and the stream function 

d 

— < u a (x 1 ,t)...u! a (x n ,t)A a (x n+1 ,t)...A a (x n+m ,t) > 

n 

^2 < ^a(xi,t)...(-e a pd a ip a (xi, t)dpAil> a (xi, t) + e oll3 d a J a (x i , t)dpA a {x u t) 
i=i 

+ vAuj a (xi,t))...uj a (x n , t)A a (x n+ i,t)...A t)> + 

n+m 

^2 < v a (xi,t)..uj a (x n ,t)A a (x n+1 ,t)..(-v aa (x i ,t)d a A a (x i ,t) + rjJ a (xi,t))..A a (x n+m ,t) > 

i=n+l 

(7) 

At each time, the probability law for both the magnetic and velocity fields is modified. We 
shall suppose that this probability law at time t after a transient period of 0(r) becomes 
a statistical field theory described by a Boltzmann weight exp(— S tjT ) 

dV{A a ^ a ) = -dAdiPeMStA A a^a)) (8) 

where the partition function Z is a normalising factor. The action St, T is a function of |> 
We shall be interested in the subsequent evolution of this probability law. 

The equations (7) can be shown to be invariant under renormalisation transformations. 
Let us denote by the scaling dimension of A and d^ the scaling dimension of if) obtained 
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from the field theory (8) . These scaling dimensions are well-defined when the field theory 
(8) is in the basin of attraction of a fixed point of the renormalisation group. In that case, 
the dimensions d,A and are the conformal dimension of the field A and ip at the fixed 
point. Let us analyse the solutions of (7). The evolution of solutions of (7) is given by 

A a (\x,\ T t) = \- d *A x -i a (x,t) 
M\x,\ T t)=\- d ^ x -i a (x,t) 

and for the action 

e -s xTttT (A a (\ x ,\ T t),j a (\x,\ T t)) = f [ d A][dil;}e~ St ^ {Ax ~ la{X ' t),i ' x ~ laiX ' t)) (10) 

where S\T tT is the action at time X T t for the solutions defined by (9). The right-hand 
side of (10) denotes a renormalisation group transformation where one integrates over all 
fluctuations of A and tp in the range [f ,a\. This renormalisation group transformation is 
carried out in momentum space where modes in the range [a -1 , Aa _1 ] are integrated out. 
We see that the probability law of the stream function remains a statistical field theory 
given by the renormalised Boltzmann weight (10). 

Let us show that (9) and (10) give the evolution of the probability law of the gauge 
field and the stream function. Let us first write (5) at time t after transforming a — > A -1 a 
and r — > X~ T r. This reads 



9tU\-i a + eapdalpX-iadpAlpX-ia = f A T 2 Aw A -i a + e al 3daJA-iadl3A\-i a 

d t A x -i a + v x -i atCt d a A x -i a = t]X t ~ 2 J x -i c 



(11) 



These equations are valid at time t in any correlation functions of the gauge field and the 
stream function. Using (9), one can substitute X dA A a and X d ^ip a for A x ~i a and ip\-i a - In 
any correlation functions, the functional integral over modes in [a -1 , Aa _1 ] can be carried 
out and involves only the Boltzmann weight. The integration of the Boltzmann weight 
gives (10). One then obtain the equations 

d X T t oj a + X 2+d ^- T e af3 d Xxa ij a d Xxfj d 2 x ^a = X 2+2dA - d o- T e af3 d Xxa J a d Xxfj A a + vd\ x u a 
d x rtA a + X 2+d ^ T v aa d Xxa A a = rjJ a . 

(12) 

These equalities are satisfied in any correlation functions at time X T t and coordinates Xx. 
They coincide with (5) provided the time rescaling exponent and the dimensions of the 
fields satisfy 

I = T 2 < i3 > 

These conditions guarantee that the evolution of the solutions of (7) is given by (9) and 
(10). The time evolution of solutions is then well-understood when the initial probability 
law of the gauge field and the stream function (8) is in the basin of attraction of a fixed 
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point of the renormalisation group. In that case, the probability law of the gauge field and 
the stream function converges to the probability law defined by the fixed point in the long 
time regime. Moreover the fixed point is a conformal field theory ^ . 

Let us analyse the physical consequences of (13). The 2-point correlation functions of 
the magnetic field and the velocity field is characterised by the scaling dimension (conformal 
dimension) cIa = d^. It is easily seen that in the long time regime, the energy spectra 
converge to 

E mag (k) ~ k 2d - +1 
£ kil #)~^ +1 

in the range k <C a -1 . The equality between the scaling dimensions of the gauge field 
and the stream function implies that these two spectra are proportional. This result 
characterises the extreme correlation between the velocity field and the stream function. 
In the following, we shall further analyse the possible fixed points and determine examples 
of energy spectra. 

In the long time regime, solutions describe a turbulent state if the viscosity and the 
resistivity become negligible. This implies that the dimension of ip has to be negative in 
order to ensure that the rescaled viscosity (respectively resistivity) u(X) = \ d ^v (respec- 
tively 77(A) = A ^77) deduced from (11) converge to zero. The fixed points are therefore 
non-unitary conformal field theories. The exponent T being positive, one deduces 

-2 < < 0. (15) 

The second inequality implies that the dimensions of the electric current J and the vorticity 
lo are positive. This guarantees that the current J and the vorticity u decrease sufficiently 
rapidly at infinity. 

Solutions can asymptotically reach a time independent non-unitary fixed point if and 
only if all the non-linear terms in (11) converge to zero. The non-linear terms can be 
readily evaluated when the cut-off is rescaled to zero using properties of short distance 
expansions. Following Polyakov, we shall suppose that the gauge field and the stream 
function become primary fields at the fixed point with conformal dimensions dA- Then the 
short distance expansion of these fields can be summarised using conformal families ([A] 
is the conformal family of A, i.e. A and all its descendent fields) 

[A][A] = [A 2 ] + ... 

MW = W + ... (16) 

MM = W + -, 

where dots mean higher order terms. Therefore the non-linear terms read 

e af3 d a ^ x -i a df3A^ x - la ~ (^-^(L^P, - L-aLiiM-ia.a + ... 
eaed a Jx-iadpA x -i a ~ (-^) dA2 ~ 2dA (L-2L 2 _ 1 - L_ 2 L* ^-i^ + ... (17) 
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where L n = / dz z n+1 T(z) are the Virasoro generators deduced from the energy mo- 
mentum tensor T of the fixed point. Using the equality between the dimensions of the 
stream function and the gauge field, these products vanish asymptotically if 

d^ 2 > 2d^ 

d A2 > 2d A (18) 
d x > 2d^p — 2. 

Any non-unitary conformal field theory satisfying these criteria is a time-independent fixed 
point of the magnetohydro dynamics equations 1 . 

In order to get more information on the fixed points, we shall be interested in the fluxes 
of the energy E = \ J d 2 x(v 2 + B 2 ) and the magnetic enstrophies M n = J d 2 xA™. In the 
Kolmogorov picture^ of turbulent steady states, these fluxes are constant. In particular, 
the constraint used by Ferretti and Yang concerns the magnetic enstrophy (n = 2). In the 
decaying case, one can expect these fluxes to vanish in the long time regime as no forcing 
is imposed. The fluxes associated to the magnetic enstrophies and the energy are ^ 1 ' 2 ^ 



RnM = ~ n J d ^ x °q * (v aa d a A a A n a 1 (x)) 

£a(q) = - J d 2 X0 q * (eapdalpadpAt/jalpa), 



(19) 



where q is an approximate delta function q (x) = J k>q d 2 k e 2nlk - x . The fluxes R n and e 
represent the transfer rate of magnetic enstrophies and energy in momentum space. In the 
asymptotic regime, the q dependence of R n , a disappears. Defining the product of n copies 
of A with itself, the product of n copies of ip and the product of n copies of A with ip ^ 

[A]...[A] = [A n ) + ... 

[#.#] = [Vn] + - (20) 
mA]...[A] = [ Xn+1 } + ... 

where dots denote higher order terms, the fluxes behave in the long time regime as 

Rn,a = a 2+ 'W-<» +1 > d * I d 2 x(L-2L 2 _ 1 — L-2L 2 _ 1 )Xn+l,X- 1 a 

, J (21) 
e a = a ^ 3 -3^ J d 2 x{L_ 2 L 2 _ x - L_ 2 L 2 _ x )^ x -r a . 

The appearance of a descendent field of Xn+i on the right hand side of (21) stems from the 
properties of the field Xn+i under parity (similarly for ^3). Notice that r n>a is a scalar and 
Xn+i is a pseudo-scalar as ip is a pseudo-scalar. The right hand side of (21) is invariant 
under parity as it is constructed using a combination of Virasoro generators which picks up 



1 The -2 was erroneously missing in Ref. [3] ® 
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a minus sign under a parity transformation. This combination is the one we had already 
encountered when calculating the non linear terms of the equations (5). The result (21) is 
valid provided 

d^) 3 — Clip I ^l/>2 ' 

In the decaying case, fluxes given by (21) have vanishing expectation values. Moreover, 
the fluxes vanish at the fixed points 

R n = 

(23) 

e = 

when the dimensions of (L_ 2 L?_ 1 — L_ 2 L?_ 1 )x n _|_i and (L_ 2 L?_ 1 — L_ 2 L^_ 1 )^3 are different 
from two. This is expected as there is no forcing on large scales. 

In general the constraints (22) are difficult to analyse. We shall construct a particular 
family of non-trivial fixed points of 2d freely decaying magnetohydrodynamic turbulence. 
The simplest examples of non-unitary conformal fields can be drawn from minimal models 
[°>. These models possess a finite number of primary fields. In order to distinguish scalars 
from pseudo-scalars, it is convenient to introduce a parity factor equal to one for scalars 
and -1 for pseudoscalars. Let us suppose that the gauge field A becomes equal to the dual 
of the stream function at the fixed point, i.e. the scalar field A = A <g> 1 differs from the 
pseudo-scalar ip = A ® — 1 only by the parity factor. This is compatible with the equality 
cIa = d-tp. Let us consider now the (p,p + 2) non-unitary minimal models (p > 1 and odd). 
One can identify A with the unique field of negative dimension such that [A] [A] = [/] + ... . 
In particular, one gets 

d A = -. (23) 

This implies that A 2 = /, ip2 = I an d x = / where I is the dual of the identity /. 
The conditions (18) are automatically satisfied. Let us now check the requirements (22). 
Using the orthogonality between primary fields, one can deduce that [A] [A] [A] = [A] . This 
entails that A n = I for n even and A n = A for n odd. Similarly one obtains ips = ip, 
Xn+i = I for n odd, and Xn+i = fp for n even. This implies that (22) is identically 
satisfied. This proves that (p, p + 2) minimal models provide a family of fixed points of 2d 
magnetohydrodynamics. It is worth noting that in this case the magnetic and the kinetic 
energy spectra are equal 

E(k) ~ fc"Mi^2) +1 . (25) 

This spectrum is always integrable. The (p, p + 2) minimal models describe finite energy 
density solutions of 2d freely decaying magnetohydrodynamic turbulence in the long time 
regime. 

Let us comment on the similarities between the decaying case and steady states of 2d 
magnetohydrodynamic turbulence. In both situations, conformal field theories are used to 
describe the probability law of the gauge field and the stream function. For steady states, 
conformal invariance is an assumption used in order to describe the scale invariant regime 
in the inertial range. In the freely decaying case, one can control the time evolution of 
solutions and prove that conformal field theories appear in the long time regime. In both 
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problems, the vanishing of (17) follows from scale invariance in a regime where the viscosity 
and the resistivity are negligible. Another similarity is the appearance of non-unitary 
theories. For steady states, this is a consequence of the constant magnetic enstrophy 
flux condition. In the free decaying case, the non-unitarity comes from the requirement 
that the renormalised viscosity v(X) = \ d ^v and the renormalised resistivity 77(A) = r\\ d * 
become negligible in the long time regime. In both cases, the fact that the viscosity is non- 
zero plays a fundamental role. This reinforces Polyakov's idea^ that non-unitary theories 
should describe turbulence as turbulence is a flux state where dissipation takes place on 
small scales and not an equilibrium state. 

Finally, there are some relevant differences between steady states and the freely de- 
caying case. First of all, solutions of 2d freely decaying magnetohydrodynamic turbulence 
do not require parity violation whereas it is easier to understand steady states of 2d mag- 
netohydrodynamic turbulence if parity is violated. We have also found that > —2 at 
a fixed point. This condition ensures that the vorticity and the electric current decay at 
infinity. Obviously, this is intimately linked to the absence of boundaries. It is certainly 
a different situation for steady states where boundary effects play a significant role. The 
same condition can be written T > 0, it is then connected to the breaking of time reversal 
invariance. Indeed, this constraint implies that solutions at later times X T t are obtained 
after averaging over small scale features at time t. Therefore in the process of evolution, 
small scale details of the gauge field and the stream function are forgotten. This implies 
that one cannot deduce solutions at time t from solutions at time X T t. In that sense, the 
behaviour of solutions of magnetohydrodynamics is irreversible. Eventually, the equipar- 
tition of magnetic and kinetic energies is not automatically satisfied for steady states. We 
have seen that the kinetic and the magnetic energy spectra are proportional in the asymp- 
totic regime of the decaying case. Unfortunately, our approach does not allow to determine 
the proportionality constant relating the kinetic and magnetic energies. 

In conclusion, we have studied the long time evolution of 2d freely decaying magneto- 
hydrodynamic turbulence assuming that the initial probability law of the gauge field and 
the stream function can be described by a statistical field theory in the basin of attraction 
of a renormalisation group fixed point. We have shown that the time evolution of this 
probability law is given by renormalisation trajectories. In that case, the long time regime 
is determined by non-unitary fixed points of the renormalistion group. In the asymp- 
totic regime, the kinetic and the magnetic energy spectra are proportional. We have then 
constructed a family of fixed points using non-unitary minimal models of conformal field 
theories. The existence of a large number of fixed points provides an explicit example of 
non-universality. The convergence towards a particular fixed point is indeed dependent on 
the initial conditions. It would be extremely interesting to understand the dynamics of 2d 
magnetohydrodynamic turbulence when the initial conditions are arbitrary. In particular, 
the influence of the infinite number of fixed points on the trajectories of solutions is a 
noteworthy problem. 
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